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Some of the most remarkable tilings and discrete quasiperiodic sets used in 
quasicrystal physics can be obtained by using strip projection method in a 
■ superspace of dimension four, five or six, and the projection of a unit hypercube 

£\J , as a window of selection. We present some mathematical results which allow one 

to use this very elegant method in superspaces of dimension much higher, and 
to generate discrete quasiperiodic sets with a more complicated local structure 
by starting from the corresponding covering cluster. Hundreds of points of these 
sets can be obtained in only a few minutes by using our computer programs. 
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. 1. Introduction 

in 

ly-j ' Quasicrystals are materials with perfect long-range order, but with no three- 

dimensional translational periodicity. The discovery of these solids in the early 
■ 1980's and the challenge to describe their structure led to a great interest in 

discrete quasiperiodic sets and their coverings ([§] and references therein). 

The diffraction image of a quasicrystal often contains a set of sharp Bragg 
peaks invariant under a finite non-crystallographic group of symmetries G, called 
the symmetry group of quasicrystal (in reciprocal space). In the case of qua- 
sicrystals with no translational periodicity this group is the icosahedral group 
Y and in the case of quasicrystals periodic along one direction (two-dimensional 
quasicrystals) G is one of the dihedral groups D$ (octagonal quasicrystals), 
Dio (decagonal quasicrystals) and D\2 (dodecagonal quasicrystals). Real struc- 
ture information obtained by high resolution transmission electron microscopy 
suggests us that a quasicrystal with symmetry group G can be regarded as a 
quasiperiodic packing of copies of a well-defined G-invariant cluster C. 

In the literature on quasicrystals the term 'cluster' has several meanings 9 . 
In the present paper, by G- cluster we mean a finite union of orbits of a finite 
group G in a fixed representation. A mathematical algorithm for generating 
quasiperiodic point sets by starting from G-clusters was proposed by author in 
collaboration with Jean-Louis Verger-Gaugry several years ago 0- It is based on 
strip projection method ( 9 and references therein) and is a direct generalization 
of the algorithm used by Katz and Duneau in 8 . The model obtained in 8 
for the icosahedral quasicrystals starts from the one-shell F-cluster C formed by 
the vertices of a regular icosahedron. The physical space is embedded into the 
superspace R 6 such that the orthogonal projections on the physical space of the 
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points 

(±1,0,0,0,0,0), (0, ±1,0, 0,0,0), (0,0, 0,0, ±1,0), (0,0,0,0,0,±1) 

are the vertices of a regular icosahedron. 

In our direct generalization, we consider only G-clusters invariant under 
inversion. If our starting G-cluster C has 2k points then we embed the physical 
space into the superspace R fe in such a way that C is the orthogonal projection 
on the physical space of the subset 

{(±1,0,0,...,0), (0,±1,0,0,...,0), (0,0,...,0,±1,0), (0,0,...,0,±1)} 

of M. k containing 2k points. One can remark that, in the case of a two-shell or 
three-shell cluster, the dimension of the involved superspace is rather high. 

Our aim is to present some mathematical results which allow one to use 
our algorithm in the case of multi-shell clusters. We show that in the case of 
a two-dimensional (resp. three-dimensional) cluster we have to compute only 
determinants of order three (resp. four), independently of the dimension of the 
superspace we use. This remark and a simple description of the window (which, 
generally, is a polyhedron with hundreds or thousands faces) have allowed us to 
obtain some very efficient computer programs for our algorithm [7] . 

In the case of a three-shell F-cluster formed by the vertices of a regular 
icosahedron, a regular dodecahedron and an icosidodecahedron we use a 31- 
dimensional superspace, the window is a polyhedron lying in a 28-dimensional 
subspace bounded by 31465 pairs of parallel faces, but we obtain 400-500 points 
in less than 10 minutes In the case of a two-shell Dio-cluster we use a 10- 
dimensional superspace, the window is bounded by 120 pairs of parallel faces, 
and we obtain 700-800 points in only one minute 



2. G-clusters 

Consider a finite group G and a fixed R- irreducible representation of G in K™. 
In the case of dihedral groups Z?2m = (a, b | a 2m = b 2 = (ab) 2 — e) we can use the 
two-dimensional representation defined by a, b : M 2 — > IR 2 

a(a, 13) = (a cos £ - (3 sin a sin £ + /3 cos 

b(a,0) = (a, -(3) [) 

and in the case of the icosahedral group Y = (a, b | a 5 = b 2 = (ab) 3 = e) the 
three-dimensional representation generated by the rotations a, b : R 3 — ► M 3 

a(a,A 7 ) = (^«-i/3±H ia+i/?+^7, 4« + ^ + il) (2) 
b(a, /3, 7) = (-a,-p, 7). 

where r = (1 + v5)/2. The set G(ai, a%, a n ) = { g(ai, ct2, a n ) I .9 £ G } 
is called the orbit of G generated by (ai, 02, a n ). If a G (0, 00) then the orbit 

Y(a, ar, 0) is formed by the vertices of a regular icosahedron, 
Y(a, a, a) is formed by the vertices of a regular dodecahedron, 
y(a,0,0) is formed by the vertices of an icosidodecahedron. 
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For each (a,/3) ^ (0,0), the orbit D 2m (a, /3) is formed by the vertices of a 
regular polygon with 2m sides. A G- cluster is a finite union of orbits of G. For 
example, C = D 2m (ai, Pi) U £>2m («2, #2) is two-shell I?2m-cluster, and 

C = Y(a, err, 0) U Y(J3, (3, (3) U Y (7, 0, 0)} 

is a three-shell icosahedral cluster. 



3. Discrete quasiperiodic sets defined by G-clusters 



Let C = {vi, i> 2 , -Ui, -«2, -v k }, where Uj = (vij,v 2 j, ...,v nj ), be a 

fixed G-cluster symmetric with respect to the origin. One can prove [21 El that 
the vectors w\, 11)2, w n , where u>i — (vn, i>i 2 , Vik), are orthogonal and have 
the same norm (which we denote by k) 



(Wi,Wj) = Vu Vjl + V i2 Vj2 + ••• + VikV-jk 



k 2 if z = j 



(3) 



if 

We identify the space K™ containing C ('physical space') with the subspace 

E = { a\Wi + a 2 w 2 + ■■■ + a„w n \ ati, a 2 , a„ el}. (4) 
of the super space R fe , and consider the orthogonal complement of E 

E 1 - = { x e R fc I (x, y) = for all y £ E }. (5) 
The orthogonal projectors corresponding to E and E arc 

E nx = (x,m)m + (x, 2*)m + ... + (x, 

\ ' K I K \ ' K I K \ ' K / K 



7T : 

n 1 - 



E J 



7T X = X — 7T X. 



(6) 



If we describe E by using the orthogonal basis {n 2 w\, k 2 w 2 , k 2 w n } then 
the orthogonal projector corresponding to E becomes 



T:R k — ► K n Vx= {(x 1 w 1 ),(x,w 2 ) 1 ...,(x,w n )). 
The projection W = -71 — L ($^) of the unit hypercube 



(7) 
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(xi,x 2 , ...,Xk) 



-\<^< \ for alHe {l,2,...,fc} 



is a polyhedron (called a window) in the (k — n)-subspace E -1 . Each (fc- 



1)- 



dimensional face of W is the projection of a (k— n— l)-face of the unit hypercube 
f2. Each (k — n — l)-face of is parallel to k — n — 1 vectors of the canonical 
basis {ei, e 2 , ek} of and orthogonal to n + 1 of them. For each n + 1 
distinct vectors e^, ej 2 , ei n+1 the number of (fc— n— l)-faces of f2 orthogonal 
to them is 2™ +1 , and the set 



.t = {xi,x 2 , 



,x kj 



Xi e {-1/2, 1/2} if i e {h, i 2 , i n +i} 
Xi = if % {ii, i 2 , i n +i} 
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Figure 1: Left: The strip S and the window W in the case of a ID physical space 
E embedded into a 3D superspace. Centre: A one-shell Dg-cluster C. Right: A 
fragment of the set Q defined by C. The nearest neighbours of any point q of Q 
belong to q + C, which is a copy of C with the center at point q. Therefore, C is 
a covering cluster for Q. 



contains one and only one point from each of them. 

There are fc!/[(n+l)! (fc-n-1)!] sets of 2" +1 parallel (fc-n-l)-faces of ft. 
In the case n = 2 these sets can be labelled by using the elements of 

1= {(21,22,23) e Z 3 | 1 < h < k- 2, + 1 <i 2 < fe- 1, i 2 + 1 < 23 < k } 

and in the case n = 3 the elements of 



(21,22,23,24) G ^ 4 



1 <h <k-3, 
12 + 1 < 23 < k — 1, 



i'l + 1 < «2 < fc - 2, 
23 + 1 < 24 < 



In the case n = 2, the vector defined by the formal determinant 
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(8) 



(9) 



for any x £ M. k . The vector y belongs to E , and since — ir^ei is a linear com- 
bination of Wi and u>2, it is also orthogonal to n ± e4, n e$, n e^. Therefore, 
y is orthogonal to the (k — 3)-faces of W labelled by (1,2,3). Similar results 
can be obtained for any 12, 23) G 1. 
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Figure 2: Left: A one-shell Di2-cluster and a fragment of the corresponding 
quasiperiodic set. Right: A fragment of the quasiperiodic set defined by a two- 
shell Z?io-cluster, obtained by using strip projection method in a ten-dimensional 
superspace. The starting cluster is a covering cluster, but for most of the points 
the occupation is extremely low. 



Consider the strip S corresponding to the window W (see figure 1) 

S = {iel l | n-t-x e w } 

and define for each (i\, i n +i) G 1 the number 
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A point x — (xi, X2, Xk) € R fc belongs to the strip S if and only if 
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for any (ii,i 2 , ...,*n+i) G 2". 

The set defined in terms of the strip projection method 



Q = v{S n z fe ) = {Vx I xesnz*} 



(10) 



(11) 



(12) 



(13) 



by using the strip S defined above and the hyper-lattice Z fc is a discrete set. If G 
is one of the groups occurring in quasicrystal physics then Q is a quasiperiodic 
set, and has all the properties of the sets obtained by projection [SI El- We have 



Ve.i = ((ei,wi), (ei,w 2 ), (ei,w n )) = {v u ,v 2i , ...,v ni ) 



(14) 
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for any i G {1, 2, k}, whence 



V({x±ei, x±e 2 , £±e fc }nS) 

C{Vx±vi,Vx±v 2 ...,Vx±v k } = Vx + C 

that is, the 'arithmetical' neighbours [5] of any point Vx <E Q belong to the 
translated copy Vx + C of C. Therefore, the starting cluster C can be regarded 
as a covering cluster ( 9' , page 16) of the point set Q. Some examples are shown 
in figure 1 and figure 2. A presentation of our algorithm in the particular case 
of a two-shell Dio-cluster can be found in [5J- 

4. Concluding remarks 

It is known that, generally, the discrete sets Q considered in the previous section 
can also be defined as multi-component model sets 00] by using superspaces of 
smaller dimension and root lattices, but generally we have to use a large number 
of very complicated windows 0] . 

The starting cluster C is a covering cluster for the set Q, but unfortunately, 
for most of the points of Q the occupation is extremely low. Therefore, our 
discrete quasiperiodic sets can not be used directly in the description of atomic 
positions in quasicrystals. 
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